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Abstract. Wc ascribe the interpretation of the twin kilohertz Quasi Periodic Oscillations (kHz QPOs) of X-ray 
spectra of Low Mass X-Ray Binaries (LMXBs) to MHD Alfven wave oscillations in the different mass density 
regions of the accreted matter at the preferred radius, and the upper kHz QPO frequency coincides with the 
Keplerian frequency. The proposed model concludes that the kHz QPO frequencies depend inversely on the 
preferred radius, and that theoretical relation between the upper frequency (u^) and the lower frequency (vi) 
is ui ~ 1^1, which is similar to the measured empirical relation. The separation between the twin frequencies 
decreases (increases) with increasing kHz QPO frequency if the lower kHz QPO frequency is more (less) than ~ 
400 Hz. 
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1. Introduction 

The launch of the X-ray timing satellite, Rossi X-ray 
Timing Explorer (RXTE), led to the discovery of Quasi 
Periodic Oscillations (QPOs) of LMXBs in their X-ray 
brightness, with frequencies ^ 10^^ — 10"^ Hz (see van der 
Klis 2000 for a recent review). Thereafter, much atten- 
tion has been paid to the kHz QPO mechanism of LMXB; 
however the proposed models are still far from explaining 
all detected data. The Z sources (Atoll sources), which are 
high (less) luminous neutron star (NS) LMXBs (Hasinger 
& van der Klis 1989), typically show four distinct types 
of QPOs (van der Klis 2000). At present, these are the 
normal branch oscillation (NBO) i^NBO — 5 — 20 Hz, the 
horizontal branch oscillation (HBO) j^hbo — 10 — 70 Hz, 
and the kHz QPOs 1^2(1^1) ^ 300 - 1300 Hz that typically 
occur in pairs in more than 20 sources, with upper fre- 
quency 1^2 and lower frequency vi. In 11 sources, nearly 
coherent burst oscillations i^burst — 270 — 620 Hz have 
also been detected during thermonuclear Type I X-ray 
bursts; these are considered to be the NS spin frequencies 
Vs or their first overtone (see, e.g., Strohmaycr & Bildsten 
2003). Moreover, the existence of a third kHz QPO has 
also been reported in three low-luminosity sources (Jonker 
et al. 2000). All of these QPOs but the burst oscillations 
have centroid frequencies that increase with the inferred 
mass accretion rate M. Furthermore, the frequencies 1^2 
and Ui, as well as the frequencies V2 and t^HBOi follow 



very similar relations in five Z sources, and the QPO fre- 
quencies of LMXBs and black hole candidates (BHC) have 
a tight and systematical correlation over three orders of 
magnitude in frequency (Psaltis et al. 1998, 1999; Belloni 
et al. 2002). 

Various theoretical models have been proposed to ac- 
count for the QPO phenomenon in X-ray binaries (for a 
review see, e.g., Psaltis 2000). In the early detection of 
RXTE, the upper kHz QPO (1^2) was simply considered to 
originate from the Keplerian orbital frequency at the pre- 
ferred radius, and the lower kHz QPO (vi) is attributed to 
the beat of this frequency with the stellar spin frequency 
i^s (Strohmayer et al. 1996; Miller et al. 1998). However, 
this beat model is inadequate, for the detected frequency 
separation (Ai^ = 1^2 — 1^1) decreased systematically with 
instantaneous M (see, e.g., van der Klis 2000). Later on, 
general relativistic effects were invoked to account for kHz 
QPOs (Stella & Vietri 1999; Stella et al. 1999; Psaltis & 
Norman 2000), which can satisfactorily explain the vari- 
ation in kHz QPO separation Ai/. Moreover, the theory 
of epicyclic parametric resonance in relativistic accretion 
disks was proposed (Abramowicz et al. 2003), where the 
twin kHz QPOs occur at the frequency of meridional oscil- 
lation and the radial epicyclic frequency in the same orbit, 
which can explain the frequency ratio 3:2 detected in black 
hole candidates. Although many other feasible ideas have 
also proposed, such as the disk seismic model (Wagoner 
1999), a two-oscillator model (Osherovich & Titarchuk) 
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has yet explained satisfactorily all observed QPO phenom- 
ena of LMXBs until now. 

In this paper, the MHD Alfven wave oscillation model 

is proposed, and its predictions and comparisons with the 
well detected sample sources are shown in the figures. 

2. The Model 

The model's idea is traced to the analogies with mag- 
netic waves in coronal oscillations in solar physics, where 
the MHD turbulence driven by Alfven wave oscillations 
occurs in the magnetic loops of the Sun's coronal atmo- 
sphere (Roberts 2000). While the dynamical details of the 
mechanisms responsible for the kHz QPOs in LMXBs are 
still uncertain, it is convenient to imagine that the MHD 
Alfven wave oscillations occur at a certain preferred ra- 
dius, where a MHD tube loop may be formed to conduct 
the accreted matter to the polar cap of star. Nevertheless, 
it is assumed that this preferred radius is a critical or 
a transitional radius where the spherical accretion matter 
with low mass density is transferred into polar channel ac- 
cretion with high mass density that follows the loop and 
accretes onto NS polar cap. Possibly, this critical transi- 
tion may give rise to MHD turbulence so that much more 
energy is liberated than at other positions. 

As a phenomenological prescription, we associate the 
twin kHz QPO frequencies with the Alfven wave oscilla- 
tion frequencies (AWOFs) at a preferred radius described 
in the Appendix, where the AWOF with the spherical ac- 
cretion mass density coincides with the Keplerian orbital 
frequency uk and the AWOF with the polar accretion 
mass density corresponds to a lower frequency. As a pre- 
liminary investigation, we are not concerned with the ac- 
tual mechanisms of producing the kHz QPOs in the X- 
ray fluxes of LMXBs (see, e.g.. Miller et al. 1998; Psaltis 
2000). Rather, our main purpose is to stress the consis- 
tence between the model and the measured kHz QPO 
data, and leave the arguments concerning the possibility 
of the mechanism to a later paper. 

At the preferred radius r, the Alfven velocity is, 
VA{f) = "^§^7 where B{r) — Bs{R/r)^ and are the 
dipole magnetic field strengths at radius r and at the sur- 
face of star with radius R, respectively. Therefore, the 
AWOF i/A is 



VA{r) _ B,{R/rf Svff 
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where the mass density of the accreted matter, p = 
M/[Svff] (Shapiro & Teukolsky 1983) is applied with 
the free fall velocity Vff = y/2GM/r = Cy/Rg/r, where 
Rs = 2GM/c^ is the Schwarzschild radius and can be ex- 
pressed as Rs — 3TO(km) = 0.3m (10km) with m = 
the NS mass in units of solar mass, and the area S rep- 
resenting the spherical area Sr or the polar cap area Sp, 
respectively. 



Sp = 47ri?2(l - cosOc) , sin^Oc = R/r = X , 



(3) 



where 6c is the open angle of the last field line to close at 
radius r. As an approximation, the polar cap area is usu- 
ally written as Sp = if i? < r (Shapiro & Teukolsky 
1983). For simplicity, it is convenient to write the two ar- 
eas by means of the scaled radius parameter X = R/r, so 
we have, 



AirR'^X-^ , Sp = 4-kR'^{1 - VT^) . 



(4) 



It is assumed that 1^2 and i^i arc from the MHD Alfven 
wave oscillations with the different accreted material mass 
densities, corresponding to the different areas Sr and 
Sp respectively, and 1/2 coincides with the Keplerian fre- 
quency at the preferred radius (see Appendix), therefore. 



1^2 = 



GM 

4^2 j,3 



= VA{Sr) = 1850(Hz)ylX3/2 , 



(5) 



with the parameter A defined as ^ = (m/i?g)^/^ and Rq = 
R/lO^cm. The physical meaning of A is clear in that 
represents the NS averaged mass density. So, by means of 
Eq.(l) with the correlation pa oc VS, we obtain the lower 
kHz QPO frequency. 



1^1 = MSp) = y2\ISp/Sr = V2X\ll-s/T^. 



(6) 



The twin kHz QPO frequencies only depend on two vari- 
ables, the parameter A and the scaled radius X, so these 
two variables are implied if the twin kHz QPO frequencies 
are detected simultaneously. Therefore, the ratio of the 
twin frequencies is obtained easily, by setting y{X) = ^, 



V2 



= y-\X)=X\^\-^Y^ 



(7) 



which is independent of the parameter A and is only re- 
lated to the parameter X. Furthermore, the twin frequency 
separation is given as follows. 



t^V = V2-Vx= V2{\ - 1 - Vl - X) . 



(8) 



The comparisons of the model's conclusions to the four 
detected sample sources, Sco X-1 (van der Klis et al. 1997), 
4U1608-52 (Mcndcz ct al. 1998a,b), 4U1735-44 (Ford et 
al. 1998) and 4U1728-34 (Mendez & van der Klis 1999), 
are shown in the figures. The relations V2 versus v\ and Ai/ 
versus vx are plotted in Fig.l and Fig. 2, respectively, and 
they show that the agreement between the model and the 
observed QPO data is quite good for the selected ranges of 
NS parameters A=0.6, 0.7 and 0.8. In Fig. 2, we find that 
Azv increases with increasing V\ if v\ < 408 (A/0.7) Hz 
and Av decreases with increasing z/i if i^i > 408 (A/0.7) 
Hz. The theoretical relation between the twin frequencies 
is derived from Eq.(5) and Eq.(6), 
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where V2k = i/2/(lkHz). If V2 < 1850A(Hz) = 
1295(g^)(Hz), then we obtain the approximated theoret- 
ical relation between the twin kHz QPO frequencies, 



A 



ui = (382M)(Hz)i/2\ = 546(Hz)(— )-i(j/2fc) 



(10) 



A similar V1-V2 empirical correlation has also been found 
for the measured kHz QPO sources (see, e.g., Stella et al. 
1999; Psaltis et al. 1998; Psaltis 2000; Psaltis and Norman 
2000). In Fig. 3, we plot the v^jvx versus v\ diagram, and it 
is found that the ratio v-i / Vi decreases with increasing i/i . 
The averaged frequency ratio for the four sample sources 
is about < y >— lA — 7 : 5, corresponding to the aver- 
aged X, < X >= 0.88. In Fig. 4, we plot 1^2/1^1 versus X, 
and find that the theoretical curve is independent of the 
parameter A, which reflects the pure geometrical scaling 
relation between the twin frequency ratio and the param- 
eter X=R/r. The X distributions for the four examples 
arc very similar (from X=G.8 to X=0.93), which implies 
that the dynamic mechanism that accounts for kHz QPO 
is intimately related to the scaled radius X and has no 
direct relation with the other physical quantities. 

In conclusion, the consistence between the model and 
the detections is robust, and the main results are sum- 
maried as follows: (1) the twin QPOs arc inversely (pro- 
portionally) related to the radii (the accretion rate); (2) 
the theoretical relation between the twin frequencies is 
vi ^ i'2- (3) the separation between the twin frequencies 
decreases (increases) with increasing the kHz QPO fre- 
quency if the lower kHz frequency is more (less) than 
400 Hz; (4) the ratio between the twin frequencies is only 
related to the scaled radius parameter X, and the homoge- 
nous kHz QPO frequency distributions of four detected 
sources indicate that these frequencies arc from regions 
close to the surface of the NS, approximately, from 1.08R 
to 1.25R. With regard to conclusion (3), for the theoret- 
ical Ai^-i^i relation with i^i < 400 Hz future detection is 
needed to confirm this prediction. 

Acknowledgements. Thanks are due to T. Belloni, M. Mendez, 
and D. Psaltis for providing the data files, and discussion with 
T.P. Li is highly appreciated. 

APPENDIX: On the preferred radius 

The Alfven radius Ra is determined by the condition 
that the magnetic ram pressure should match the plasma 

momentum energy density (Shapiro & Teukolsky 1983), 

where Ra = 1.9 x \Q'^{cm)B'^J^ M~^^'^ m-^l^ r]^^\ and 
Sgg the surface field in unit of 10^ Gauss, M17 the ac- 
cretion rate in unit of lO^'' g/s, 



^^ = \p{Ra>\Ra). 



(11) 



or equivalently, Wa(-Ra) = 



_ B(Ra) _ 



Vi;^ - v{Ra\ where va 
and V represent the Alfven velocity and the free fall ve- 



and p(r) = Pff{r) = M/[5w//(r)] with the area factor 
S = 47rr^. In other words, the Alfven velocity equals the 
free fall velocity at the Alfven radius. If the area factor S 
decreases (the polar cap area for instance) , the mass den- 
sity p will increase and the Alfven velocity will decrease. If 
we set v{r) = Vk{r) and p{r) = Pff{r), then we will obtain 
a preferred radius where the Alfven velocity matches the 
Keplerian flow velocity. 



va{R, 



'SI 



B{R. 



SI) 



^4ttp{Rsi) 



= Vk{Rsi) ■ 



(12) 



So, the radius Rsi can be obtained to be Rsi = 
22/7i?^ = 1.22Ra ~ (B«/Mi/2)4/7^ tijg Alfven 

wave oscillation frequency matches the Keplerian fre- 
quency, z^a(-Rsi) = i^k{Rsi)', we prefer to call it the quasi 
sonic-point radius as a distinction of the sonic-point radius 
by Miller et al. (1998) and Lai (1998). 
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Fig. 1. vi versus v\ plot. The horizontal axis is the lower kHz 
QPO frequency v\ and the vertical axis is the upper kHz QPO 
frequency v^. The kHz QPO data of four detected sample 
sources are plotted. The model presents a good consistence 
with the measured data for the NS averaged mass density pa- 
rameters A=0.6, 0.7 and 0.8, which are shown in the three 
theoretical curves from bottom to top. 
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Fig. 3. The ratio between the twin kHz QPO frequencies versus 
the lower kHz QPO frequency. The three theoretical curves rep- 
resent the NS mass density parameter conditions with A=0.6, 
0.7 and 0.8 from bottom to top. The data of four sample sources 
are plotted. 
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Fig. 2. Ai^ versus v\ plot. The horizontal axis is the lower 
kHz QPO frequency v\ and the vertical axis is the twin kHz 
QPO frequency separation Aj^. The kHz QPO data of four 
detected sample sources are plotted. The model presents a good 
consistence with the measured data for the NS averaged mass 
density parameters A=0.6, 0.7 and 0.8, which are shown in the 
three theoretical curves from bottom to top. 
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Fig. 4. The ratio between the twin kHz QPO frequencies ver- 
sus the parameter X = R/r. The horizontal axis is X and the 
vertical axis is the ratio j=V2/i'i- It is found that the theoret- 
ical curve is independent of the NS mass density parameter A. 
The two horizontal lines represent y=5:4 and y=5:3, from bot- 
tom to top, respectively, which cover the ratios of the detected 
kHz QPO data of the four samples. 
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Abstract. We ascribe the interpretation of the twin kilohertz Quasi Periodic Oscillations (kHz QPOs) of X-ray 
spectra of Low Mass X-Ray Binaries (LMXBs) to MHD Alfven wave oscillations in the different mass density 
regions of the accreted matter at the preferred radius, and the upper kHz QPO frequency coincides with the 
Keplerian frequency. The proposed model concludes that the kHz QPO frequencies depend inversely on the 
preferred radius, and that theoretical relation between the upper frequency {1^2) and the lower frequency (vi) 
is fi ~ 1^1, which is similar to the measured empirical relation. The separation between the twin frequencies 
decreases (increases) with increasing kHz QPO frequency if the lower kHz QPO frequency is more (less) than ~ 
400 Hz. 
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1. Introduction 

The launch of the X-ray timing satellite, Rossi X-ray 
Timing Explorer (RXTE), led to the discovery of Quasi 
Periodic Oscillations (QPOs) of LMXBs in their X-ray 
brightness, with frequencies ~ 10~^ — 10^ Hz (see van der 
Klis 2000 for a recent review). Thereafter, much atten- 
tion has been paid to the kHz QPO mechanism of LMXB; 
however the proposed models are still far from explaining 
ah detected data. The Z sources ( Atoh sources) , which are 
high (less) luminous neutron star (NS) LMXBs (Hasinger 
& van der Klis 1989), typically show four distinct types 
of QPOs (van der Klis 2000). At present, these are the 
normal branch oscillation (NBO) i^bo — 5 — 20 Hz, the 
horizontal branch oscillation (HBO) vubo — 10 — 70 Hz, 
and the kHz QPOs 1^2(^1) ^ 300 - 1300 Hz that typically 
occur in pairs in more than 20 sources, with upper fre- 
quency 1^2 and lower frequency Vi. In 11 sources, nearly 
coherent burst oscillations fc^urst 

~ 270 - 620 Hz have 
also been detected during thermonuclear Type I X-ray 
bursts; these are considered to be the NS spin frequencies 
Vs or their first overtone (see, e.g., Strohmayer & Bildsten 
2003). Moreover, the existence of a third kHz QPO has 
also been reported in three low-luminosity sources (Jonkcr 
et al. 2000). All of these QPOs but the burst oscillations 
have centroid frequencies that increase with the inferred 
mass accretion rate M. Furthermore, the frequencies i>2 
and vi, as well as the frequencies 1/2 and i/ubo, follow 
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X-rays: stars 

very similar relations in five Z sources, and the QPO fre- 
quencies of LMXBs and black hole candidates (BHC) have 
a tight and systematical correlation over three orders of 
magnitude in frequency (Psaltis et al. 1998, 1999; Belloni 
et al. 2002). 

Various theoretical models have been proposed to ac- 
count for the QPO phenomenon in X-ray binaries (for a 
review see, e.g., Psaltis 2000). In the early detection of 
RXTE, the upper kHz QPO (2^2) was simply considered to 
originate from the Keplerian orbital frequency at the pre- 
ferred radius, and the lower kHz QPO (fi) is attributed to 
the beat of this frequency with the stehar spin frequency 
Vg (Strohmayer et al. 1996; Miher et al. 1998). However, 
this beat model is inadequate, for the detected frequency 
separation {Av = V2 — Vi) decreased systematically with 
instantaneous M (see, e.g., van der Klis 2000). Later on, 
general relativistic effects were invoked to account for kHz 
QPOs (Stella & Vietri 1999; Stella et al. 1999; Psaltis & 
Norman 2000), which can satisfactorily explain the vari- 
ation in kHz QPO separation Af. Moreover, the theory 
of epicyclic parametric resonance in relativistic accretion 
disks was proposed (Abramowicz et al. 2003), where the 
twin kHz QPOs occur at the frequency of meridional oscil- 
lation and the radial epicyclic frequency in the same orbit, 
which can explain the frequency ratio 3:2 detected in black 
hole candidates. Although many other feasible ideas have 
also proposed, such as the disk seismic model (Wagoner 
1999), a two-oscihator model (Osherovich & Titarchuk) 
and the photon bubble model (Klein et al. 1996), no model 
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has yet explained satisfactorily all observed QPO phenom- 
ena of LMXBs until now. 

In this paper, the MHD Alfven wave oscillation model 
is proposed, and its predictions and comparisons with the 
well detected sample sources are shown in the figures. 

2. The Model 

The model's idea is traced to the analogies with mag- 
netic waves in coronal oscillations in solar physics, where 
the MHD turbulence driven by Alfven wave oscillations 
occurs in the magnetic loops of the Sun's coronal atmo- 
sphere (Roberts 2000). While the dynamical details of the 
mechanisms responsible for the kHz QPOs in LMXBs arc 
still uncertain, it is convenient to imagine that the MHD 
Alfven wave oscillations occur at a certain preferred ra- 
dius, where a MHD tube loop may be formed to conduct 
the accreted matter to the polar cap of star. Nevertheless, 
it is assumed that this preferred radius is a critical or 
a transitional radius where the spherical accretion matter 
with low mass density is transferred into polar channel ac- 
cretion with high mass density that follows the loop and 
accretes onto NS polar cap. Possibly, this critical transi- 
tion may give rise to MHD turbulence so that much more 
energy is liberated than at other positions. 

As a phenomcnological prescription, we associate the 
twin kHz QPO frequencies with the Alfven wave oscilla- 
tion frequencies (AWOFs) at a preferred radius described 
in the Appendix, where the AWOF with the spherical ac- 
cretion mass density coincides with the Keplerian orbital 
frequency vk and the AWOF with the polar accretion 
mass density corresponds to a lower frequency. As a pre- 
liminary investigation, we are not concerned with the ac- 
tual mechanisms of producing the kHz QPOs in the X- 
ray fluxes of LMXBs (sec, e.g.. Miller et al. 1998; Psaltis 
2000). Rather, our main purpose is to stress the consis- 
tence between the model and the measured kHz QPO 
data, and leave the arguments concerning the possibility 
of the mechanism to a later paper. 

At the preferred radius r, the Alfven velocity is. 



Sp = 4itR^{1 - cos6a) , si-n?6a = R/r = X 



(3) 



where 6c is the open angle of the last field hne to close at 
radius r. As an approximation, the polar cap area is usu- 
ally written as Sp = if i? <C r (Shapiro & Teukolsky 
1983). For simplicity, it is convenient to write the two ar- 
eas by means of the scaled radius parameter X = R/r, so 
we have. 



4TrR^X-'^ , Sp = 47rJ?2(l - Vl - X) . 



(4) 



It is assumed that 1^2 and vi are from the MHD Alfven 
wave oscillations with the different accreted material mass 
densities, corresponding to the different areas Sr and 
Sp respectively, and 1/2 coincides with the Keplerian fre- 
quency at the preferred radius (see Appendix), therefore. 



2^2 



GM 



PAiSr) = 1850(Hz)AX 



3/2 



(5) 



with the parameter A defined as A = {ui/Rq)^^'^ and Rq = 
R/lO^cm. The physical meaning of A is clear in that 
represents the NS averaged mass density. So, by means 
of Eq.(??) with the correlation i/a oc \/S, we obtain the 
lower kHz QPO frequency. 



Vl = Va{Sp) = V2 



^SplSr = V2X\ll-y/T^ . 



(6) 



The twin kHz QPO frequencies only depend on two vari- 
ables, the parameter A and the scaled radius X, so these 
two variables are imphed if the twin kHz QPO frequencies 
are detected simultaneously. Therefore, the ratio of the 
twin frequencies is obtained easily, by setting y{X) = 



1^2 



y-'{X)=X^Jl-VT^.. 



(7) 



which is independent of the parameter A and is only re- 
lated to the parameter X. Furthermore, the twin frequency 
separation is given as follows. 



Mil 



/47rp 



, where B{r) = B,{Rlrf and B, are the = ^2 - vi = i'2{l - X\Jl 



VA{r) 

dipole magnetic field strengths at radius r and at the sur- 
face of star with radius R, respectively. Therefore, the 
AWOF VA is 



X) 



(8) 



VAir) 



VA{r) _ B^jR/rY Svff 



27rr 



27rr 



A-kM 



\fs, 



(1) 



where the mass density of the accreted matter, p = 
M/[Svff] (Shapiro & Teukolsky 1983) is apphed with 
the free fall velocity Vff = ^JlGMjr = c^jR^/r, where 
Rs = 2GM /c^ is the Scliwarzscliild radius and can be ex- 
pressed as Rg ~ 3m(km) = 0.3TO(10km) with m = 
the NS mass in units of solar mass, and the area S rep- 
resenting the spherical area Sr or the polar cap area Sp, 
respectively. 



The comparisons of the model's conclusions to the four 
detected sample sources, Sco X-1 (van der Klis et al. 1997), 
4U1608-52 (Mendez et al. 1998a,b), 4U1735-44 (Ford et 
al. 1998) and 4U1728-34 (Mcndcz & van der Klis 1999), 
are shown in the figures. The relations V2 versus vi and Av 
versus v\ are plotted in Fig.l and Fig. 2, respectively, and 
they show that the agreement between the model and the 
observed QPO data is quite good for the selected ranges of 
NS parameters A=0.6, 0.7 and 0.8. In Fig. 2, we find that 
A/y increases with increasing vi if vi < 408 (A/0.7) Hz 
and Af decreases with increasing vi if vi > 408 (A/0.7) 
Hz. The theoretical relation between the twin frequencies 
is derived from Eq.(??) and Eq.(??), 



47rr^ 



(2) 



vi = 629iHz)A 




1.&5A' 



(9) 
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where V2k 



z/2/(lkHz). If V2 < 1850A(Hz) 



1295((j^)(Hz), then we obtain the approximated theoret- 
ical relation between the twin kHz QPO frequencies, 



1^1 = (382M)(Hz)z.2\ 



546(Hz)(— )-i(z.2.)^ (10) 



A similar 1^1-1^2 empirical correlation has also been found 
for the measured kHz QPO sources (see, e.g., Stella et al. 
1999; Psaltis et al. 1998; Psaltis 2000; Psaltis and Norman 
2000). In Fig. 3, we plot the 1^2/1^1 versus diagram, and it 
is found that the ratio 1^2/1^1 decreases with increasing Vi. 
The averaged frequency ratio for the four sample sources 
is about < y >= 1.4 = 7 : 5, corresponding to the aver- 
aged X, < X >= 0.88. In Fig. 4, we plot 1^2/ 1^1 versus X, 
and find that the theoretical curve is independent of the 
parameter A, which reflects the pure geometrical scaling 
relation between the twin frequency ratio and the param- 
eter X=R/r. The X distributions for the four examples 
are very similar (from X=0.8 to X=0.93), which implies 
that the dynamic mechanism that accounts for kHz QPO 
is intimately related to the scaled radius X and has no 
direct relation with the other physical quantities. 

In conclusion, the consistence between the model and 
the detections is robust, and the main results are sum- 
maried as follows: (1) the twin QPOs are inversely (pro- 
portionally) related to the radii (the accretion rate); (2) 
the theoretical relation between the twin frequencies is 
ui ~ 1^1; (3) the separation between the twin frequencies 
decreases (increases) with increasing the kHz QPO fre- 
quency if the lower kHz frequency is more (less) than ~ 
400 Hz; (4) the ratio between the twin frequencies is only 
related to the scaled radius parameter X, and the homoge- 
nous kHz QPO frequency distributions of four detected 
sources indicate that these frequencies are from regions 
close to the surface of the NS, approximately, from 1.08R 
to 1.25R. With regard to conclusion (3), for the theoret- 
ical Au-ui relation with < 400 Hz future detection is 
needed to confirm this prediction. 

Acknowledgements. Thanks are due to T. Belloni, M. Mendez, 
and D. Psaltis for providing the data files, and discussion with 
T.P. Li is highly appreciated. 

APPENDIX: On the preferred radius 

The Alfven radius Ra is determined by the condition 
that the magnetic ram pressure should match the plasma 
momentum energy density (Shapiro & Teukolsky 1983), 
where Ra = 1.9 x 10'^' {cm)D'll^ M~^^^''m-^/'^Rl^^\ and 
Bs8 the surface field in unit of 10*^ Gauss, Mn the ac- 
cretion rate in unit of lO^'^ g/s, 



^^ = \piRAyiRA)^ 



(11) 



or equivalently, va(Ra) = f^^^^ = v(Ra), where va 

and V represent the Alfven velocity and the free fall ve- 
locity, v{r) =Vff{r) = ^IGMjr = c^Rsjr = V2vkir) 



and p{r) = Pff{r) = M /[Svff{r)] with the area factor 
S = Anr^. In other words, the Alfven velocity equals the 
free fall velocity at the Alfven radius. If the area factor S 
decreases (the polar cap area for instance) , the mass den- 
sity p will increase and the Alfven velocity will decrease. If 
we set v{r) = Vk{r) and p(r) = Pff{r), then we will obtain 
a preferred radius where the Alfven velocity matches the 
Keplerian fiow velocity. 



va{R. 



SI) 



B{Rsi 



^yinpiRsi) 



Vk{R. 



SI) 



(12) 



So, the radius Rgi can be obtained to be Rsi = 
2'^I'^Ra = 1.22Ra ~ (Bs/MV2)4/7^ where the Alfven 
wave oscillation frequency matches the Keplerian fre- 
quency, va{Rsi) = i^k{Rsi)] we prefer to call it the quasi 
sonic-point radius as a distinction of the sonic-point radius 
by Miller et al. (1998) and Lai (1998). 
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Fig. 1. vi versus v\ plot. The horizontal axis is the lower kHz 
QPO frequency v\ and the vertical axis is the upper kHz QPO 
frequency v-x. The kHz QPO data of four detected sample 
sources are plotted. The model presents a good consistence 
with the measured data for the NS averaged mass density pa- 
rameters A=0.6, 0.7 and 0.8, which are shown in the three 
theoretical curves from bottom to top. 
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Fig. 3. The ratio between the twin kHz QPO frequencies versus 
the lower kHz QPO frequency. The three theoretical curves rep- 
resent the NS mass density parameter conditions with A=0.6, 
0.7 and 0.8 from bottom to top. The data of four sample sources 
are plotted. 
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Fig. 2. Ai/ versus v\ plot. The horizontal axis is the lower 
kHz QPO frequency v\ and the vertical axis is the twin kHz 
QPO frequency separation Af. The kHz QPO data of four 
detected sample sources are plotted. The model presents a good 
consistence with the measured data for the NS averaged mass 
density parameters A=0.6, 0.7 and 0.8, which are shown in the 
three theoretical curves from bottom to top. 
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Fig. 4. The ratio between the twin kHz QPO frequencies ver- 
sus the parameter X = R/r. The horizontal axis is X and the 
vertical axis is the ratio j=i'2 /i^i ■ It is found that the theoret- 
ical curve is independent of the NS mass density parameter A. 
The two horizontal lines represent y=5:4 and y=5:3, from bot- 
tom to top, respectively, which cover the ratios of the detected 
kHz QPO data of the four samples. 
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